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Starting with the Vlasov-Bohzmann equation for a binary fluid mixture, we derive an equation 
for the velocity field u when the system is segregated into two phases (at low temperatures) with 
a sharp interface between them, u satisfies the incompressible Navier-Stokes equations together 
' with a jump boundary condition for the pressure across the interface which, in turn, moves with 

, a velocity given by the normal component of u. Numerical simulations of the Vlasov-Boltzmann 

■ equations for shear flows parallel and perpendicular to the interface in a phase segregated mixture 

support this analysis. We expect similar behavior in real fluid mixtures. 



^ . PACS numbers: 64.75.+g, 47.70.Nd, 68.05.-n 
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When a fluid is quenched from a high temperature to a low temperature inside the miscibility gap it evolves from 
I 1 its initial homogeneous state which is now unstable into a final equilibrium state consisting of two coexisting phases 
[ separated by an interface. The phase segregation process, involving hydrodynamical velocity fields u(x, t), is more 
complicated than the corresponding diffusive processes in quenched binary alloys In a seminal work Siggia 
pointed out that in the late stages of phase segregation sharp interfaces develop between domains of the two phases 
, and the coarsening process is mainly due to hydrodynamic flows driven by these elastic interfaces. Following Siggia 
^ ■ there has been much work to investigate the time evolution of this coarsening process P, 0, Isj , and this is still going 

^: on[|g. 

' O ■ There has been much less work on the derivation of a quantitative description of the velocity field involving the 
5 ' interface between two phases of macroscopic size, a problem of considerable interest in its own right 1^. To do this 
one should start from a microscopic description of the transition region in which the composition varies smoothly on 
the microscopic scale and end up with boundary conditions on the "membrane" representing this transition layer on 
the macroscopic scale (the sharp interface limit). This is a difficult task for the case in which the two phases have 
' widely different densities, viscosities, etc., since a reasonable microscopic description valid for both phases is hard to 
find. 

To overcome this difficulty we consider here a situation in which both the equilibrium and kinetics are as simple as 
possible. This is a binary fiuid model introduced in |^ which provides a description of phase segregation kinetics at the 
00 mesoscopic level via the Vlasov-Boltzmann (VB) equations. These equations can be obtained from the microscopic 
dynamics (in a convincing but not rigorous manner) and in turn allow a derivation, in a suitable limit, of bulk 
, hydrodynamical equations [^ . The hydrodynamic equations we obtain for the bulk mixture are similar but not 
identical to the heuristic macroscopic equations previously proposed for describing phase segregation in an energy 
and momentum conserving system [ pT| . In the present paper we consider the case in which the ffuid is already 
well segregated with sharp interfaces between the different phases. We derive hydrodynamic equations involving the 
G interface motion, directly from the VB equations and carry out direct simulations for the time evolution of simple 
hydrodynamical fiows, in the presence of such an interface. The resulting free boundary problem was considered 
phenomenologically, in the context of interface oscillations, by Harrison p2| . We compare our results, with good 
O ' agreement, with the behavior based on the linearized version of the hydrodynamic equations, recently considered 
O for studying coarsening dynamics |^ . Both the derivation and the numerical analysis are new and we believe that the 
^ , results have, as is the case of the Navier-Stokes equations of a one component fluid, a universal validity even though 
their rigorous derivation from microscopic models can only be carried out currently from the Boltzmann equation 
. In fact preliminary calculation using the Enskog equation instead of the Boltzmann equation (more appropriate 
. for a dense fluid) leads to the same hydrodynamics. 
■ " ' The model we consider is a mixture of two species of hard spheres with diameter d and unit mass labeled by 
1,2, with a long range positive potential £~^V{£~^r) between particles of different species; £ is the range of the 
interaction. The latter causes phase segregation at temperatures lower than a critical temperature Tc, into phases I 
and II: I is rich in species 1 and II in species 2 Let n be the total particle density. It was shown in that, 
when n(P << 1, while £ >> d, this system is well described by coupled VB kinetic equations for the one-particle 
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distributions fi{r,v,T),i = 1,2 

drh + V ■ Vrh + F, ■ V,/,; = J[/,; , /i + /a] (1) 

where Fi{r,T) = —WrVi{r,T) with Vi the self-consistent Vlasov potentials Vi(r,r) — J^g^-^ £^^V{£^^\r — 
r'\)nj{r', T)dr' , i ^ j, nj{r, t) — Jj^., /,(r, v,T)dv, and J[/, g] is the Bohzmann colhsion operator for hard spheres [^ . 
We are interested in an equihbrium state at a fixed total density in which the system is split into two regions, one con- 
sisting of phase 1, and the other of phase 11. These regions are separated by a transition layer whose local form can be 
obtained by considering particular 1-d stationary 'solitonic' (actually kink) solutions of set fi = pi{z)M{v) where 
M{v) is the Maxwellian with unit mass, zero mean velocity and temperature T, Af(v) = (27rfcBT)~^/^ exjp[—v'^/2kBT]. 
Moreover Pi{z) — (l)i(z/t) with 4>i(ci) smooth functions satisfying the equations 

ftsriog0,(q) -f / dq'V{\q - q'\)c^,{q') = (2) 

where V{q) — J^2 dwV{y^lp~\-~w^) and d is a constant. For T < Tc Eqs. (§) have non-homogeneous solutions with 
asymptotic values and pj^ as q ^ ±oo, corresponding to the densities in phases 1 and II respectively: by symmetry 
Pi — p¥ J pi = p¥ ■ The solutions are not explicit but their existence and properties are studied in |^ and |]l5j. These 
are similar to those observed in realistic mixtures, such as oil and water. 

We are interested in studying hydrodynamical flows when the width of the interface, which is of order £, is small 
compared to the mean free path A which is in turn small compared to L, the characteristic length of the domains 
occupied by the pure phases. Setting r ~ e^^x,T ~ f^^t, where e = X/L ^ £/X and x, t are the macroscopic 
position and time, we wish to study the small e behavior of a solution ffjx, v, t) — fi{e~^x, v, e~'^t) of Eqs. (1), in the 
incompressible regime, corresponding to low hydrodynamic velocity |14|] . We assume far from the interface a solution 
of the form 

OO 

f^{x, V, t) = -p,{x)M{v) + ^ ^gl(x, V, t) (3) 

s=l 

where pi{x) is either pi or pl^ depending on which region x is in. Note that by symmetry Pi + pi = pV + p¥ — P: 
i.e. the total density is the same in both phases but varies in the interface region. Putting = /i + /|, we have for 
the hydrodynamic velocity field 

p^W vj%x,v,t)dv ^ e u{x,t) + 0{e^), 

SO that the Mach number is of order e (incompressible regime). Far from the interface the temperature T and the 
densities are constants, coinciding for T <Tc with the asymptotic values of the solitonic solutions. 

To understand what happens near the interface consider the situation in which the interface is flat. Then /^(r, i;, r) = 
Pi{z)M{v) with pi{z) = (j)i{z/£) a stationary solution of the Vlasov-Boltzmann equations. If the interface is not flat 
this is not exactly true, but because of the tendency of the solitonic profile to keep its form and just translate, one 
expects that the solution is locally close to the solitonic profiles multiplied by the Maxwellian. On the other hand, 
the force due to the surface tension acts on the particles of the fluid which start to move with some average velocity 
u and the surface moves accordingly, so that its normal velocity is at any point u ■ TV, N being the normal to the 
surface. To implement this picture we write the solution near the interface as 

oo 

f^{z,x,v,t) ^ p,{z)M{v) + J2^"9t{z,x,v,t). (4) 

s=l 

The notation is as follows: given a point x we call z — e^^d{x,rt), where d(x,rt) is the distance of x from the 
interface Ft, and x the component of x tangential to the interface; Pi{z) — (j)i(z/£) is the solitonic solution centered 
on z = 0. Since we are so close to the interface we can assume that locally it looks almost flat, so that the solution 
near the interface has a weak dependence on x. Therefore the profile interpolating between the values of the densities 
in the bulk on the two sides of the interface should be well approximated by the 1-d solitonic profile Pi{z), the lowest 
order term in Eq. (^. 

To obtain a solution of Eqs. (1) we have to put these expressions in the equations (after space-time rescaling) and 
match the two expansions. This can be done in a consistent way and it is possible to compute the terms in the series 
(the long computation will be reported elsewhere) , but the question of the convergence of the series is open (we expect 
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they are asymptotic). The result is that in the hmit e the velocity field u{x,t) is divergence free and solves the 
incompressible Navier-Stokes (INS) equation 

df U + {u ■ V)u + = i^Au (5) 

with the kinematic viscosity u obtained from the Boltzmann equation as in Moreover, u is continuous across the 
interface Tt whose normal velocity is given by 

vr,{x) =u{x,t) ■ N{x,t) (6) 

while the pressure is discontinuous at the surface and satisfies Laplace's law 

ip+-p-)=aK, (7) 

Here (p_) p+ is the pressure on the side of Ff (not) containing the normal N; K is the mean curvature of Tf and a 
is the surface tension given in terms of solitonic profiles as 

E '^V{<l^l')U<l')dqdq' (8) 

The equation for m is a complicated initial boundary value problem with free boundary Tt ■ The solutions of INS 
in the two domains separated by Tt are coupled by the jump condition for the pressure. Since the solution is not 
explicit, we compare our numerical results with the linearized version of this problem, studied in JlSf . 

To obtain numerical solutions of the VB equations we use the method introduced in |^ to simulate the VB dynamics 
at the particle level, based on coupling the direct simulation Monte Carlo algorithm for the short-range interaction with 
the grid-weighting method for the long-range repulsion. The results were obtained using about 5 x 10^ particles, in a 
cube (LxLxL) with periodic boundary conditions. The unit of length is the mean free path A = (2^ Trpd^)"^, and the 
unit of time is the mean-free time between collisions r = A/c, c = {2kBT)i . For the repulsive inter-species potential 
we use, as in [||, V{q) = a7r~2 exp(— g^) (note that the mean-field critical temperature is given by ksTc = pa/2), 
with a range of interaction £ — 0.4A. All quenches were performed at total particle density p, with pd^ « 0.01. We 
focus on the behavior of a simple initial fiow perturbation, Uz{0) = uocos(fcy), with k = 2Tmk/L, Uk = 1,2,3,..., 
outside and inside the coexistence region; for all simulations we set uq — 0.1c. For a homogeneous system this shear 
wave perturbation leads to a time dependent flow Uz{t) — u{t) cos{ky), where u{t) — u{0) exp(— z^/c^t), ly — rj/ p and rj 
is the viscosity. 

We first study the behavior of the above initial perturbation in the homogeneous region of the phase diagram at 
T/Tc = 1.5. The initial velocity profile decays exponentially with very high accuracy and we extract a kinematic 
viscosity i/ which agrees with the Boltzmann gas viscosity |l^. We then turn to the study of the system in the 
presence of interfaces, at two temperatures, Ti/Tc = 0.51 and T2/TC = 0.33. For each of these temperatures the initial 
configuration contains two interfaces (due to periodic boundary conditions), situated in the {xy} plane a distance 
L/2 apart, that separate domains of the two equilibrium phases for the particular temperature. The static structure 
of these interfaces has been described in The thickness of the interface is of order £. 

For an initial flow perturbation that is parallel with the interfaces, u(0) = Ux{0)A = uq cos(fcy)x, where x is the unit 
vector in the x-direction, the behavior of the velocity profile is virtually identical with the one in the homogeneous 
region. The situation is however very different if the initial perturbation is perpendicular to the interfaces, i.e. u(0) = 
Mz(0)z — uocos{ky)z. While Uz{t) is still very well represented by u{t)cos{ky), u{t) no longer decays exponentially. 
To fully characterize this situation we also look at the position of the interface. The order parameter ip — (pi — P2)/ P 
is calculated at time as a function of z and y. We find that ip{z,y,t) is very well fitted by the profile (^otanh[(z — 
Zo{y,t))/£] 0, and identify zo{y,t) with the position of the interface. Furthermore, zo{y,t) is well represented by 
Zo{y,t) = A(t)cos{Ky). The velocity amplitude u{t) and the interface amplitude A{t) for the temperature Ti and 
Ufe = 1 are shown in Fig. 1, along with fits with the functional forms expected for an overdamped harmonic oscillator: 
A{t) = Mo exp(— 7t) sinh(a't)/a;, u{t) = exp(— 7i)[— 7sinh(ci;t)/w -I- cosh(ci;t)]. 

We find that these forms can reasonably reproduce the results up to the highest simulated, Uk — 9. In Fig. 2 
we show the amplitudes of the interface oscillations A{t) up to — 3, along with the corresponding best fits, for 
temperature Ti(top) and r2 (bottom). We note that the Uk = 1 behavior at T2 is qualitatively different from the other 
cases, showing damped oscillations. The functional forms that describe this result are: A{t) — uq exp{—^t) sm{ujt) / ui , 
u(t) = uoexp(— 7t)[— 7sin(a;i)/ijj -f cos(a;t)]. 

These simulation results can be compared with those obtained from the linear perturbation analysis of a thin, 
initially fiat interface with surface tension a separating fiuids of matching density p and viscosity 77 = vp. This type of 
analysis has been carried out in ||l3] in the full space. The idea is to look for flow solutions u^^y.zit) oc e'^* cos{kxX+kyy) 
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following a small velocity perturbation of wave number k. This yields a Uz decaying exponentially as a function of the 
distance from the interface. We carried out a similar analysis matching the periodic boundary and initial conditions 
in which Uz is independent of z away from the interfaces. Under these conditions and taking into account the 
discontinuity of the normal stresses due to surface tension, we can write, following jl^ (see e.g. Chapter X, Eq. 28) 

pdtuz = v^u, + ^sHd^ + d^)6zs]Siz - z,) (9) 

where Zs is the position of the interface s, Szs is the displacement from the initial position and dSzs/dt = Uzs is 
the velocity of the interface. However, we regard the above equation integrated over the torus, along z, as the more 
fundamental one if Uz is z independent. This yields 

pdtUz = v^uz + L-^^s<y[dl + dl)5zs (10) 

Using Uz = Uzs and assuming Uz{t) oc e^* cos(fcy), which is consistent with the initial conditions, we get p? p + ^ik^rj + 
L~^2ak'^ — 0. This equation is characteristic of a damped oscillator. Introducing the length scale Lg = pja and 
time scale Tg = (?■ ja^ ^ this becomes 

M*+Ai*fc*+-^*"^2/c2 = (11) 

where /z* = /xT,, /c* = kL^ and = L/Lg. Using k = 27rnfc/L we find that the above equation predicts damped 
oscillations for Uk < n-l = (2L*) 2 /n and overdamped behavior (no oscillations) otherwise. For our binary fluid model 
the surface tension can be related to the interface profile (see Eq. 8) and easily estimated jT^ , and therefore we can 
calculate n^. We find that for Ti/Tc = 0.55, n° = 0.83, so that for all the oscillations should be overdamped, 
in agreement with the simulation results. For T2/TC — 0.33, n° = 1.39, so the = 1 case should exhibit damped 
oscillations, while the others should be overdamped; this is again in agreement with the simulations. We also calculate 
the coefficients 7 and uj that describe the behavior of the velocity profile and compare with the simulation results, see 
Fig. 3. We find good agreement for small k (nk) and increasing deviations as k increases; the relative behavior of 7 
and to also appears to be qualitatively different from the simulations for large k. This is as expected, since for larger fc, 
i.e. smaller wavelengths, we are far from the hydrodynamic regime. Our hydrodynamic picture is valid if kX = 27rnfee, 
£/\ and u/c are of order £ = j;, which is in our simulations 0.028. Better agreement is found in simulations with 
£/X — 0.2 (and e = .056), see Fig. 3, which indicates that the accuracy of the equations we derive increases when the 
interface becomes sharper. 

The work S.B. was performed under the auspices of the U. S. Department of Energy by University of California 
Lawrence Livermore National Laboratory under Contract No. W-7405-Eng-48. The work of J. L. was supported by 
NSF Grant DMR 98-13268, AFOSR Grant AF 49620-01-1-0154, DIMACS and its supporting agencies, the NSF under 
contract STC-91-19999 and the N.J. Commission on Science and Technology, and NATO Grant PST.CLG.976552. 
R.M. partially supported by MIUR, Cofin-2000 and INFM. 



[1] For a review see, e.g., J. S. Langer, in Solids far from equilibrium, edited by C. Godreche (Cambridge University Press, 

New York, 1991). 
[2] E.D. Siggia, Phys. Rev. A 20, 595 (1979). 
[3] W.R. White, P. Wiltzius, Phys. Rev. Lett. 75, 3012 (1995). 
[4] F. Perrot et al., Phys. Rev. E 59, 3079 (1999). 
[5] V.M. Kendon et al., Phys. Rev. Lett. 83, 576 (1999). 

[6] F.J. Soils, M. Olvera de la Cruz, Phys. Rev. Lett. 84, 3350 (2000); Phys. Rev. Lett. 85, 4408 (2000). 
[7] H. Furukawa, Phys. Rev. Lett. 85, 4407 (2000). 

[8] See, e.g., D.M. Anderson, G.B. McFadden, A.A. Wheeler, Annu. Rev. Fluid Mech. 30, 139 (1998). 

[9] S. Bastea, J.L. Lebowitz, Phys. Rev. Lett. 78, 3499 (1997). 
[10] S. Bastea, R. Esposito, J.L. Lebowitz, R. Marra, J. Stat. Phys., 101, 1087 (2000). 
[11] K. Kawasaki, T. Koga, T. Kawakatsu: Mat. Res. Soc. Symp. Proc. 237, 87 (1992). 
[12] W.J. Harrison, Proc. London Math. Soc. 6, 396 (1908). 

[13] S. Chandrasekar, Hydrodynamic and Hydromagnetic Stability (Clarendon Press, Oxford, 1961), Chap. X. 

[14] A. De Masi, R. Esposito, J. L. Lebowitz, Comm. Pure and Appl. Math. 42, 1189 (1989). 

[15] E. Carlen, M.C. Carvalho, R. Esposito, J. L. Lebowitz, R. Marra, in preparation. 

[16] C. Cercignani, The Boltzmann equation and its applications ( Springer- Verlag, New York, 1988). 

[17] S. Bastea, S. Puri, J.L. Lebowitz, Phys. Rev. E 63, 041513 (2001). 



5 




T ' 1 ' r 




I I I I I I I I I 

200 400 600 800 



FIG. 2: Interface amplitudes A{t) for Ti/Tc = 0.51 (top) and Tg/Tc = 0.33 (bottom), Uk = 1,2,3 (top to bottom) (dotted 
lines) and corresponding fits (see text) (solid lines). 



7 




FIG. 3: Coefficients 7 and a;(inset): from simulations (open symbols) at Ti with I = 0.4A (circles) and t = 0.2A (triangles); 
oscillator equation (see text) (filled symbols). 



